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We calculate the Debye and Meissner masses of photons and gluons for a spin-one color supercon- 
ductor in the polar and color-spin-locked phases as well as for a spin-zero color superconductor in 
the 2SC and color-flavor-locked phases. A general derivation for photon-gluon mixing is provided in 
terms of the QCD partition function. We show qualitatively and quantitatively which of the gauge 
bosons attain a mass via the Anderson-Higgs mechanism. Contrary to the spin-zero phases, both 
spin-one phases exhibit an electromagnetic Meissner effect. 

PACS numbers: 12.38.Mh,24.85.+p 



I. INTRODUCTION 

Sufficiently cold and dense quark matter is a color superconductor. Due to asymptotic freedom, one-gluon exchange 
is the dominant interaction in a dense quark system. The one-gluon exchange interaction is attractive in the color- 
antitriplet channel. Similar to electron Cooper pairing in ordinary superconductors [1], this leads to the formation 
of Cooper pairs at the quark Fermi surface [2, 3] (for a review, see [4]). Cold and dense matter could occur in the 
interior of compact stellar objects, such as neutron stars [5]. 

In ordinary superconductivity, because of spontaneous symmetry breaking of the electromagnetic gauge group 
U(l) em , the photon attains a mass via the Anderson-Higgs mechanism. Physically, this mass corresponds to screening 
of the magnetic field in superconducting matter. One of the consequences is the electromagnetic Meissner effect, i.e., 
the expulsion of magnetic fields. In order to investigate the Anderson-Higgs mechanism in a color superconductor, 
one has to take into account the color gauge group SU(3) C in addition to the electromagnetic U(l) em . The symmetry 
breaking pattern of SU(3) C x U(l) em determines which gauge bosons (gluons and photons) attain a mass. The order 
parameter for quark Cooper pairing defines the specific symmetry breaking pattern. In this paper we provide a general 
treatment of the influence of color superconductivity on color and electromagnetic fields and apply it to four different 
phases of superconducting quark matter. 

First, we consider a system of two quark flavors. Here, two quarks can form Cooper pairs in the color-antitriplet, 
flavor-singlet channel [2], commonly called the 2SC phase. The Cooper pairs in this phase carry total spin zero, J = 0. 
Moreover, in the 2SC phase the quarks of one color are supposed to remain unpaired. 

Second, we study three quark flavors in the color-flavor-locked (CFL) phase [6]. In this phase, the condensate 
consists of pairs in the color-antitriplet, flavor-antitriplet channel. All colors and flavors participate in pairing and 
the condensate is invariant under joint color-flavor rotations. Also in this phase the Cooper pairs carry spin zero. 

In the third and fourth phases considered here, quarks form spin-one Cooper pairs, J = 1 [2, 7, 8]. For symmetry 
reasons, a total spin J = 1 is required for Cooper pairs of quarks which carry the same flavor. The simplest case is 
a system of only one quark flavor. Another possibility is a many-flavor system where each flavor separately forms 
Cooper pairs. In the polar phase [9, 10], the spin of the pair points in a fixed spatial direction and, as in the 2SC 
phase, one color remains unpaired. In the color-spin-locked (CSL) phase [2], each spatial direction is associated with 
a direction in color space. The case of spin-one Cooper pairing of quarks is in some respect similar to atomic pairing 
in superfluid Helium 3 [11]. Therefore, similar phases occur in superconducting quark matter, e.g., the CSL phase 
corresponds to the "B phase" in Helium 3. 

Since the order parameter has a particular color-flavor-spin structure in each phase, the pattern of symmetry 
breaking differs from phase to phase, which consequently have to be investigated separately in order to identify the 
gauge bosons which attain a mass through the Anderson-Higgs mechanism. In some phases, namely the 2SC, CFL, and 
polar phases, it turns out that the physical fields are not the original gluon and photon fields but linear combinations 
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of these, i.e., there is mixing between photons and gluons [12, 13]. The mixing angle is the analogue of the Weinberg 
angle in the standard model of electroweak interactions. 

Our paper is organized as follows. In Sec. II we discuss the QCD partition function for massless quarks in the 
presence of gluon and photon gauge fields. We include a diquark condensation term to account for color supercon- 
ductivity. The goal of this section is to determine the gauge boson propagators. This fundamental study provides the 
basis for the following sections. It is general in the sense that we consider gauge fields for arbitrary four-momenta P 
while afterwards we focus on the limiting case P — > in order to discuss screening of static and homogeneous electric 
and magnetic fields. 

In Sec. Ill, we use a group-theoretical argument, namely the invariance of the gap matrix under some subgroup 
of SU(3) C x U(l) em , to show which of the gauge bosons attain a mass via the Anderson-Higgs mechanism. In this 
section we determine in a qualitative manner if electric and magnetic gauge fields are screened. Moreover, we derive 
the mixing angles describing the rotation of the gauge fields. 

In the following two sections we extend the results of Sec. Ill by performing a quantitative calculation of the Debye 
and Meissner masses for all gauge bosons. These masses correspond to the inverse screening lengths for electric and 
magnetic fields, respectively. They are obtained from the longitudinal and transverse components of the polarization 
tensors in the zero-energy, low-momentum limit. The first of these two sections, Sec. IV, deals with the technical 
details of the calculation. In its first part, Sec. IV A, we perform all calculations that do not depend on the special 
color-superconducting phase. A part of these calculations, namely the integrals over the absolute value of the quark 
momentum, are deferred to App. A. In the second part, Sees. IV B - IV E, we specify the results for the above 
mentioned four color-superconducting phases, i.e., we derive the particular expressions for the relevant components of 
the polarization tensors. Readers who are not interested in the technical details may skip this section. All results and 
a discussion are given in Sec. V. We present all screening masses for photons and gluons which yield the mixing angles 
and the masses for the (physically relevant) mixed, or rotated, gauge bosons. One of the points we discuss in this 
section is in which cases these mixing angles differ from the ones obtained in Sec. III. We conclude and summarize 
our study in Sec. VI. 

Our convention for the metric tensor is g^ v — diag{l, —1, — 1, —1}. Our units are h = c = ks = 1. Four- vectors are 
denoted by capital letters, K = = (ko, k), and k = |k|, while k = k/fc. We work in the imaginary-time formalism, 
i.e., T/VJ2k = TJ2 n I rf 3 k/(27r) 3 , where n labels the Matsubara frequencies u) n = ik . For bosons, u n = 2mrT, for 
fermions, uo n = (2n + l)nT. 



In this section we provide the theoretical basis for the calculation of the Debye and Meissner masses in a color 
superconductor. We clarify how "mixing" of gluon and photon fields has to be understood. The results of this section 
are the final form for gluon and photon propagators, Eq. (38), and the mixing angle between gluon and photon fields, 



We start from the grand partition function of QCD for massless quarks in the presence of gluon fields A% and a 
photon field A^, 



II. MIXING OF GLUONS AND PHOTONS 



Eq. (43). 



A. Gluon and photon propagators 




(1) 



Here, Sa is the action for gluon and photon fields. It consists of three parts, 



Sa = S F 2 + S g f + Sfpg , 



(2) 



where S g f and Sfpg are the gauge fixing and ghost terms, respectively, and 




(3) 



is the gauge field part. The space-time integration is defined as J x = J^ T dr J v d 3 x, where T is the temperature and 
V the volume of the system. The field strength tensors F£ v — d^A^ — d u A^ + gf ahc A h ^A c v correspond to the gluon 
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fields A", while = d^A v — d v A^ corresponds to the photon field A^. The functional Z q [A] is the grand partition 
function of massless quarks in the presence of gluon and photon fields and a chemical potential fi. It is given by 



Z q [A] = Jv^V^ exp J $ (i -fd^ + M7o + 9 -fA^T a + e j^A^Q) ,/> 



(4) 



where T a are the generators in the fundamental representation of the gauge group of strong interactions, SU (3) c , and Q 
is the generator of the electromagnetic U (l) em . The coupling constants for the strong interaction and electromagnetism 
are denoted by g and e, respectively. The quark fields ip are spinors in color, flavor, and Dirac space. 

In order to take into account Cooper pairing of the quarks we include a diquark source term. This has been done 
for instance in Refs. [2, 14]. We generalize this treatment by also taking into account the photon field. Then, after 
integrating out the fermion fields, the partition function is given by [14] 



/ 



VAexp 



Sa + IttMs-' + a^) 



(5) 



where the trace runs over space-time, Nambu-Gor'kov, color, flavor, and Dirac indices. The sum over a now runs 
from 1 to 9, where A 9 ^ = A^ is the photon field. We also defined the 2x2 Nambu-Gor'kov matrices 

f diag( 5 ^T a , -g 7 "T a T ) for a = 1, . . . , 8 , 
r^diagfeO^ (6) 
V ' { diag(e7^Q,-e 7 ^Q) for a = 9. 

S = S(X,Y) is the quasiparticle propagator in Nambu-Gor'kov space. In Eq. (5), we did not explicitly keep the 
fluctuations of the order parameter field, cf. Eq. (26) of Rcf. [14]. This is possible, since we are only interested in the 
gluon and photon propagator. Nevertheless, it should be kept in mind that, like in any theory with spontaneously 
broken gauge symmetry, these fluctuations mix with the (4-)longitudinal (unphysical) components of the gauge field. 
As usual, using a suitable choice of 't Hooft gauge, the fluctuations can be decoupled from the gauge field. Then, in 
unitary gauge, one finds that the gauge boson propagator is explicitly (4-)transverse [14], in accordance with general 
principles [15, 16]. 

Expanding the logarithm in Eq. (5) to second order in the gauge fields, we obtain 

^Tr InGS- 1 + Aft*) ~ ^Tr lnS" 1 + \ Tr[A« f £ S\ ±Tr[A^ f £ S A h v t\ S\ . (7) 

Following Ref. [14], the sum of all terms which are quadratic in the gauge fields will be denoted by Si. S2 does not 
only contain the pure gluon and photon terms but also two terms which mix gluon and photon fields. In order to 
perform the trace over space-time, we introduce the Fourier transforms 

S(X,Y) = Lj2^ K<X - Y) S(K) , (8a) 

K 

A%{X) = J2e- ip - x A;(P), (8b) 
p 

where we used translational invariance, S(X,Y) = S(X — Y). Then we obtain 

Tr[A;(X) f £ S(X, Y) Al (Y) f I S(Y, X)} 

K,P 

= -^E A £(- p ) n -6( p )^( p ) ' ( 9 ) 

p 

where the trace now runs over Nambu-Gor'kov, color, flavor, and Dirac indices and where we defined the polarization 
tensor 

Kl(P)^\^Y. T ^ S ^ K ) t bS{K - P)] . (10) 

K 



X.Y 
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The Nambu-Gor'kov quasiparticle propagator in momentum space is denned as 

(") 

with the quasiparticle propagators G ± (K) and the "anomalous" propagators E^ 1 (K) , 

G ± (K) = {[G±}- 1 (K)-<f>*(K)G*(K)$ ± (K)}- 1 , 5±(K) = -G* (K)$± (K)G± (K) . (12) 

Here, G^(K) = (j^K^ ± /i7o) _1 is the free propagator. In these expressions, we have put the regular self-energies (cf. 
Ref. [4]) to zero since, to the order we are computing, they do not influence the results for the polarization tensors. 
Following Ref. [10], we write the gap matrix $ + for color-superconducting quark matter, where ultrarelativistic quarks 
form Cooper pairs in the even-parity channel, as 



k • 



(13) 



e=± 



Here, (jf{K) is the gap function, A k = (1 + e7o7 • k)/2 are projectors onto states of positive or negative energy, and 
Alk is a matrix in color, flavor, and spin space that determines the color-superconducting phase. One can choose Aik 
such that [A^k, A k ] = 0. Furthermore, <I>~ (K) = 70 (<& + {K))^ 70. For the color-superconducting phases considered 
here, 7o-M k .Mk7o = M^M^. In this case, 



G ± (K) = {GT(K)]- 1 J2 E ^ A k £ 



e=± r=l,2 



1,2 
«, 



(0 e 



(14) 



where are projectors onto the eigenspaces of the hermitian matrix 7oAl k Alk7o- In all phases considered here, 
there are two eigenvalues of this matrix, denoted by Ai and A2. They appear in the quasiparticle excitation energies 



,21 1/2 



The projectors are given by 



^1.2 _ 7o- A 4- M k7o - A 2 ,i 



Al,2 — A21 



The anomalous propagators can be written as 



S+(^) = -^7oMk7onA k < 



k 2 



~k,r < 



2 ' 



E-(K) = -J2MlV r k At 



k 2 



-k,r v 



As in Ref. [14], we introduce a set of complete, orthogonal projectors for each 4- vector P M = (po,p), 



QT = tT 



OS" 



N 2 



QT 



p 2 



(15) 



(16) 



(17) 



(18) 



With = (p p 2 /P 2 ,pQp/P 2 ), the projector projects onto the one-dimensional subspace that is (4-)orthogonal to 
P 11 but (3-)parallcl to p. The operator projects onto the one-dimensional subspace parallel to P^. Consequently, 
Qi" projects onto a two-dimensional subspace that is (4-)orthogonal to both P^ and N^. Furthermore, this subspace 
is (3-)orthogonal to p. With the additional tensor 



Q1 V = N^P V + P"N" 

we can decompose the polarization tensor [16], 



(19) 



(20) 



(In the notation of Rcfs. [14, 16], Qi = A, Q 2 = B, Q 3 = E, Q 4 
coefficients W ab (P) of the projection operators are given by 



C.) Since for % = 1,2,3, Q^Qa^ = 0, the 



K b (P) = 



i = 1,2,3 



(21) 
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The remaining coefficient corresponding to <2 4 is 



,4 



n - (p) " is^r ■ (22) 

The explicit forms for W ab (P) are given in Ref. [14], Eqs. (42) and (43). Employing the decomposition of the 
polarization tensor in Eq. (9), we obtain 

s * = -^EE^(- p )2rn: h (pK(p) . (23) 



2T 
p 



Now we add the free gauge field term 



s ™ = -^E^(- p )( p V v - P"P")^(P) = -i^A«(-P)P 2 (Qr + QTWAP) • (24) 



We obtain 



2 4 



*(o ) = _iv; 

9 T ■ 

P I, i=l i=3 



5 2 + sis = — = e ^(- p ) E er + n ^( p )] + E 2rni 6 (p) 4(p) . (25) 



Since we finally want to read off the gluon and photon propagators, we have to transform this expression in two 
ways. First, concerning the Dirac structure it is necessary to get rid of the term proportional to Q4 which mixes 
the longitudinal mode (3-parallel to p) with the unphysical mode (4-parallel to P^). Then, inverting the inverse 
propagator becomes trivial, because it is just a linear combination of the complete, orthogonal projectors Qi, Q2, Q3. 
Second, in order to obtain the physical modes we have to diagonalize the resulting 9x9 matrices which, after 
eliminating Q 4 , will replace S ab P 2 + II^P), i = 1,2, and II* 6 (P), i = 3, 4 in Eq. (25). 
We first write S2 + 5^2 as 



F 2 

S2 + s$ = -~ E { E a H-p) + n «f( p )] A U p ) + A 3°" (-P)n^(P)4 M (P) 

p 1=1 

+^(-P)JV„n^(P)P I ,A^(P) + A°"(-P)P M n^(P)JV v A^ (P)} , (26) 



where A^(P) = Q^A^(P) are the gauge fields projected on the subspace corresponding to Qi. Now one can "unmix" 
the fields A%{P) and A%(P) by the following transformation of the unphysical field component A3 (P), which does 
not affect the final result since we integrate over all fields in the partition function, 

ATi-P) — A?{~P) A?(-P)N„ni b (P) [n 3 (P)]-; P^ , (27a) 
A^(P) — A^(P) - P„ [n^P)]" 1 n^(P) ^A-(P) . (27b) 

After this transformation, one is left with quadratic expressions in the projected fields. The transformation modifies 
the term corresponding to i = 2, 



S2 + s$ = -~E^(- p )(fin*-6P 2 + ni 6 (P)] 
p 

+Qr{^P 2 + n2 b (p)-p 2 ^n 4 QC (P) [n 3 (P)] c -; n^p)} + Qfn 3 ab (p))^(P) . (28) 



Before we do the diagonalization in the 9-dimcnsional gluon-photon space, we add the gauge fixing term S g f. We 
choose the following gauge with gauge parameter A, 

s g f = -^Y, A ^- p ) p,xpVA ^ p ) = -^^E^(- p ) p2 2f A "( p ) • ( 29 ) 

p p 

This looks like a covariant gauge but, including the fluctuations of the order parameter, which we did not write 
explicitly, it is actually some kind of 't Hooft gauge, cf. Eq. (50) of Ref. [14]. Moreover, had we fixed the gauge prior 



6 



to the shift (27) of the gauge fields, we would have to start with an expression which is non-local and also involves 
the (3-)longitudinal components of the gauge field. 
Adding the gauge fixing term (29) to Eq. (28) leads to 



1 V 3 



with 



S 2 + S^ + S 9f = ---J2J2 Aa ^- p ^i V@ ib(P)Ai(P) , (30) 

P i=l 



6 ab P 2 +nl b (P) for i = l, 

s ab P 2 + nl b (p) P 2 N*iit(p) [n 3 (P)] -j n 4 h (p) for i = 2 , (31) 

S ab kP 2 + Kb( P ) for i = 3. 



In order to obtain the physical modes, we have to diagonalize the 9x9 matrices 6* ft (P). Since S^ h (P) is real and 
symmetric, diagonalization is achieved via an orthogonal transformation with a 9 x 9 matrix Oi(P), 

S2+Sf 2 +S af =-\^Y.ilKA- P )Q"^aa{P)K^P) > (32) 

P i=l 

where 

A^(P)=0?(P)A^P) (33) 



are the rotated gauge fields and 



©L(p)-of(p)0L(p)or(p) (34) 



are diagonal matrices. The index i in A® ^ has a different origin than in Af^ introduced in Eq. (26). For A^, it 
indicates that for each i = 1,2,3 one has to perform a separate diagonalization. For Af^ it denoted the projection 
corresponding to the projector Q? v ' . 

Note that the orthogonal matrix Oi{P) depends on P M . For energies and momenta much larger than the super- 
conducting gap parameter the polarization tensor is explicitly (4-)transverse, II 3 = 0, and diagonal in a, b. In this 
case, Oi(P) — > 1. Consequently, the gauge fields are not rotated. However, in the limit po = 0, p — > it is known 
that gluons and photons mix at least in the 2SC and CFL phases [12], Oi(P) ^ 1. Thus, the mixing angle between 
gluons and photons, which will be discussed in Sec. II B, is in general a function of P^ and interpolates between a 
nonvanishing value at p = 0, p — > and zero when po,p — > oo. Note also that the orthogonal matrix Oi(P) depends 
on i, i.e., it may be different for longitudinal and transverse modes. We comment on this in more detail below. 

From Eq. (32), we can immediately read off the inverse propagator for gluons and photons. It is 



A- i :^)=EsreL^)- (35) 
i=i 

From the definition of Q l aa (P) we conclude 

P 2 + IT aa (P) for 2 = 1,2, 
&UP) = { (36) 
ip 2 +HL(P) for 2 = 3, 



where 



Of(P)W bc {P)Or(P) for i = l,3, 



Ka(P)={ , , , (37) 

o?{P) |n 2 c (p) - p 2 n 2 uI(p) [n3(p)] d ; nj c (P)} ot{p) for i = 2 . 

In the case po = 0, using Eqs. (42) and (43) of Ref. [14] one realizes that the extra term involving II 3 and II 4 for i = 2 
vanishes. Thus, in this case one only has to diagonalize the original polarization tensors n^ b . 
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Finally, we end up with the gauge boson propagator 

Ka (P) = = QT + = Q2 " + -~ Qf • (38) 

P2 +n i a (F) 1 P 2 +II 2 Q (P) 2 p2 + An 3 a(P) 3 

Setting the gauge parameter A = 0, we are left with the transverse and (3-)longitudinal modes, in accordance 
with general principles [15, 16]. The static color and electromagnetic properties of the color superconductor arc 
characterized by the Debye masses TfiD.a and the Meissner masses ifiM.a which are defined as 

fh^ a = - hm ITL(0,p) = - lim n°° a (0,p) , (39a) 
p— >o p— >o 

m 2 ^ = hm nL(0, p) = \ lim(^ - j^n^O, p) . (39b) 

p^O 2 p^O 

Since the orthogonal matrices Oi(P) are regular in the limit po = 0, p —> 0, the masses can also be obtained by first 
computing lim p ^o (0j p) ana - then diagonalizing the resulting 9x9 mass matrix. In Sec. IV we use this method 
to compute fh 2 D a and rh 2 M , since the diagonalization of the matrix H^(P) for arbitrary P M is too difficult. 



B. The mixing angle 

In this section we investigate the structure of the orthogonal matrices Oi (P) which diagonalize the gauge field part 
of the grand partition function. In general, the matrices Oi{P) mix all gluon components among themselves and 
with the photon. However, in the limit po — 0, p — > it turns out that the only non-zero off-diagonal elements of 
the tensor W ab = lim p ^o IT ab (0, p) are IPg 9 = n| = ni 8 . Physically speaking, gluons do not mix among themselves 
and only the eighth gluon mixes with the photon. In this case, Eq. (37) reduces to the diagonalization of a 2 x 2 
matrix. Consequently, the diagonalization is determined by only one parameter and the (nontrivial part of the) 
transformation operator reads 

( cose one A m 

1 — sin fc>, cos Vi 1 



The new fields are 



A 8 ^ = cosdiAl + sm6 i A li , (41a) 
A^i = - sin 6i Al + cos 9, A^ . (41b) 



The new polarization functions (eigenvalues of 11^) are 



ft* 8 = n| 8 cos 2 0i + 2n 87 sin t cos 6, + n; 7 sin 2 t , (42a) 
n 77 = II 88 sin 2 9i - 2n 87 sin 6i cos 6 t + cos 2 6 t . (42b) 

The mixing angle Oi is given by 

211* 

tan2(9 4 = — ^_ . (43) 

TP — TP 
11 88 11 77 

If [II g7 ] 2 = II 88 II 77 , the determinant of II* 6 is zero, which means that there is a vanishing eigenvalue. In this case, 
we have 

cos 2 0. = . . , (44) 

n 88 + n 77 

and the new polarization tensors, Eqs. (42), have the simple form 

n 88 = n 88 + n 77 , n; 7 = o. (45) 



Physically, the vanishing polarization tensor for the rotated photon corresponds to the absence of the Meissner effect 
for i = 1, or the absence of Debye screening for i = 2. 
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III. SYMMETRIES AND MIXING 



In Section II B we have discussed the mixing angle 9i that combines the electromagnetic field and the eighth gluon. 
The index i refers to the projectors Q^ v and thus, in general, the mixing angle could be different for transverse modes, 
z = l, and longitudinal modes, i = 2. 

In this section, we present a method to determine the mixing angles via a simple group-theoretical consideration. 
As we know from ordinary superconductors, the non-zero electric charge of a Cooper pair leads to the Meissner effect 
and thus corresponds to a non-vanishing Meissner mass m,M,~f for the photon. Besides this magnetic screening there 
is also electric screening of photons described by the photon Debye mass mc, 7 . Of course, in a color superconductor, 
in addition to the electric charge we also have to take into account the color charge. The group-theoretical method 
applied in the following allows us to investigate whether there exists a (new) charge which generates an unbroken 
symmetry. In this case, a Cooper pair is neutral with respect to this charge, and consequently one expects that there 
is neither a Meissner effect nor Debye screening. The new charge is a linear combination of electric and color charges. 
Correspondingly, the associated gauge field is a linear combination of the photon and the eighth gluon field, which, in 
turn, defines the mixing angle. The group-theoretical method only allows to identify a new charge, and thus does not 
distinguish between electric Debye or magnetic Meissner screening. Consequently, the mixing angles for longitudinal 
(electric) and transverse (magnetic) modes deduced by this method are identical, 0\ = 9 2 = 0. 

A Cooper pair is neutral with respect to a charge, if the gap matrix $ is invariant under a transformation generated 
by an operator Q corresponding to this charge. In general, the gap matrix is a matrix in color, flavor, and spin space, 

$ = cj) ijk e ? <g> e\ ® e{ , (46) 

where (pijk is the order parameter and <g> ® e'[ is a basis for the color£g)flavor(g>spin representation of the group 

SU(3) C x SU(Nf)f x SU{2)j. Since Q generates a subgroup of the product group of strong and electromagnetic 
interactions, SU(3) C x U(l) em , it can be found by the following invariance condition, 

$ = (fife X gemMgc X > ( 47 ) 

where g c <G SU(3) C and g em e U(l) em act on ef and ej, respectively. With g c = exp(zK a T a ), g em — exp(mQ), for 
infinitesimal transformations this condition can be written as 

= 4>i jk [«„ {T a e1 + e?T a T ) ® e{ <8> e{ + e\ ® in (q e f j + ejQ T ) ® ejj] , (48) 

and finding the charge Q is equivalent to determining the real coefficients n a , n, where a — 1,...,8. Using the 
representations and the order parameters for each phase as given in Table I we obtain 



2SC 
CFL 

CSL 
polar 



= in a (T a e c 3 + e^Tj) <g> e f <g> e J + i n e% ® (Q e f + e f Q) ® e J , (49a) 

= i n a (T a e1 + e^Tj) ® e{ <8> e J + i k e) ® (q + e^Q^j ® e J , (49b) 

= in a (T a ej + e^Tj) (gie 7 ® e] + 2inqe L j ® e f O e/ , (49c) 

= iK a (T a el + elT^)®e f ®e'l + 2inqel®e f ®ei . (49d) 



With the special form of the charge generator Q, given in the third column of Table I, Eqs. (49) yield conditions for 
the coefficients n a and n, 

2SC: Ki,K2, n 3 arbitrary, K4 = . . . = k 7 = 0, n s = -=n , (50a) 

v3 

2 

CFL: K\ = . . . = K7 = 0, k 8 = —j= K j (50b) 

v3 

CSL: ki = . . . = k 8 = k = , (50c) 
polar: «i, k 2 , «3 arbitrary, K4 = . . . = k 7 = 0, Kg = — 2V3qn . (50d) 

The fact that K\,K2, «3 are arbitrary in the 2SC and polar phases means that there is an unbroken SU(2) C symmetry 
under which the order parameter is invariant. This is obvious because in these phases Cooper pairs carry anti-blue 
color charge which is not seen by the generators of this SU(2) C , since they operate exclusively on the red-green 
subspace. The relation between n% and n determines a new U(l) symmetry, generated by 

g = Q + 77T 8 , (51) 
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representation 


order 
parameter 


Q 


generators of 
residual group 


?7 


2SC 


3 C (g> 1/ ® lj 


<5« 


diag(2/3, -1/3) 


Ti,T 2 ,T 3 ,Q + 77T 8 


-1/73 


CFL 


3 C ® 3/ ® lj 




diag(-l/3,-l/3,2/3) 


Q + r/Tg 


2/V3 


CSL 


3 C ® 3j 










polar 


3 C ® 3./ 




9 


ri,T 2 ,T3,Q + r7T 8 


-2^9 



TABLE I: Representations, order parameters <^ijfc, charge generators Q, generators of the residual symmetry group, and 
coefficients r) for all color-superconducting phases considered here. For the CSL and polar phases, the electric charge of the 
quark is denoted by q. The last two columns result from Eqs. (49). 



where r\ is given in the fifth column of Table I for the different cases considered here. In the CSL phase, the fact that 
all k's vanish means that there is no unbroken residual symmetry. 

In order to deduce the mixing angle 9 we rewrite the covariant derivative, = <9 M — igA a ^T a — ieA^Q, in terms of 
the new charge generator Q and the linearly independent generator 

T$ =T$ + \Q , (52) 

which belongs to the broken part of the group; A is a real constant to be determined below. We also have to replace 
the gauge fields and A^ by the rotated fields (41) and the associated coupling constants g, e by new coupling 
constants g, e. Thus, we demand the identity [12] 

gAlT s + eA^Q = gAffi + eA^Q . (53) 
Inserting Eqs. (41) for the rotated fields and the definitions of Q and T 8 we determine 

e 2 

g = g cos 9 , e — e cos 9 , A = — rj— , (54) 



and the mixing angle 



cos 2 9 = f . (55) 



Since j»e, the mixing angle in all three cases is small. Thus, according to Eqs. (41), the gluon almost remains a 
gluon and the photon almost remains a photon and therefore it is justified to call the rotated gauge bosons the new 
gluon and the new photon. 

In the 2SC, CFL, and polar phases, the new charge is neither Debye- nor Meissner-screened. From the argument 
presented at the end of Section II B we therefore expect that the polarization tensor for the rotated photon vanishes 
in the zero-energy, zero-momentum limit, 

n^ 7 2 - . (56) 

As we shall see in the following section, this argument is not quite correct for the 2SC and polar phases, as it neglects 
the effect of the unpaired blue quarks on the screening of electric and magnetic fields. In the CSL phase, on the other 
hand, there is no residual U(l) symmetry and we conclude that all gauge fields experience the Meissner effect and 
Debye screening, since in this case there does not exist any charge with respect to which a Cooper pair is neutral. 
Concerning color screening, one concludes from Table I that in the 2SC phase and the polar phase only the gluons 
4-8 are screened while in the CFL and CSL phases all gluons are screened. 



IV. CALCULATION OF THE POLARIZATION TENSORS 



In this section, we calculate the polarization tensors II^(P) given in Eq. (10) in the limit po = 0,p — > 0. In this 
case, the Debye and Meissner masses, cf. Eqs. (39), arc obtained from the coefficients of the first two projectors in 
the decomposition (20). They will be calculated in the next section. Here, we first derive a general expression for the 
polarization tensor that holds for all different phases. Then we insert the order parameters of the 2SC, CFL, polar, 
and CSL phases, and show the results for each phase separately. 
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A. General structure of 11^(0) 



We start from Eq. (10) and first perform the trace over Nambu-Gor'kov space, 

Kl(P) = ~Y.{^aG + {K)T v b G + {K - P)]+Tr^ a G-(K)r b G-(K - P)] 



K 



+Tr[r^ S-(K) Y v b E+(K - P)] + Tr[^ E+(K) ^ E~(K - P)]} , 



(57) 



where the traces now run over color, flavor, and Dirac space. In the following we first consider the traces with the 
quark propagators and afterwards investigate the traces containing the anomalous propagators S 1 * 1 . 

In order to find the results for the former, we first perform the Matsubara sum. This is completely analogous to 
the calculation of Ref. [17]. The only difference is our more compact notation with the help of the projectors P£, cf. 
Eq. (14). Thus, abbreviating K\ = K , K 2 = K — P , and fcj = |kj| for i = 1,2, we conclude 

T^Tr[r^G+(^ 1 )r^G+(^ 2 )] = EE^KTo^AkTr^o^A^] v™>+(k u k 2 ,p ) , (58a) 



fco 



(58b) 



where (cf. Eq. (40) of Ref. [17]) 



rs,+ (u u \ - ( ni >ri l - n 2,s) (1 - ni, r )n 2 , a \ , . 
UeieI(«l) fe 2,Po) = - ; : (1 - N hr - A 2 ) 



(1 - ni, r )(l - n 2 , s ) ni ir n 2 , s 



Po - ei, r + e 2;S Po + ei,r-e2,« 



(JVi, r - AT 2 , S ) , 



C 2 ( fc l> fc 2,Po) 



Po + Cl,r + e 2,s Po - ei,r - e 2,s / 
ni iT .n 2 , s (1 - ni )7 .)(l - n 2iS )' 



Here, we abbreviated 



e »,r — £ k,,r 



rii, 



Pq - Cl,r + £2,s 



Pa + £l,r - £2,s 



exp(ej, r /T) + 1 



(i = l,2) . 



Note that, for p = 0, we have 



iC:(*i,to,o) = «^;(fci,fe,o) = < s e2 (fci,fc 2 ,o) 



(59a) 



(59b) 



(60) 



(61) 



Next we discuss the traces containing the anomalous propagators. Again the Matsubara sum is completely analogous 
to the calculation in Ref. [17]. Therefore, using Eq. (17), we obtain 



r^Tr[r^-(x 1 )rrs+(^ 2 )] = ^^[ T " M l V ^\ T hoMu 2 7oKA 

feo ei,e 2 r,s 

T ]T Tr [rtH+^O ^ E-(K 2 )] = E E Tr [ZlaM^oV^AC^lKA: 



-e 2 
k 2 



w r e s ie2 (k u k 2 ,p ) , (62a) 
to.Po) , (62b) 



where (cf. Eq. (93) of Ref. [17]) 

w r e s ie2 {ki,k 2 ,po) 



_ <Pl,r(P2,s 



4ei jr e 2iS 



1 



Po + ei,r + £2,s Po - ei,r - £2, 

1 1 



(1 - N hr - N 2 , s ) 



Pa - £l,r + e 2 : s Po + £l,r - £2,; 



(63) 



Here, 
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is the gap function on the quasiparticle mass shell given by the excitation branch ko = £%, r - In the derivation of Eqs. 
(58) and (62), we assumed that, in the 2SC and CFL phases, the chemical potentials for all Nf quark flavors are the 
same, fii = . . . = /xjv, = At. In this case, the functions v and w only depend on the single chemical potential \x and can 
thus be factored out of the flavor trace. In the cases where quarks of the same flavor form Cooper pairs, i.e., in the 
polar and CSL phases, our formalism also allows for the treatment of a system with Nf > 1 and different chemical 
potentials fi n , n= 1, . . . , Nf. Then, v and w depend on the quark flavor through [i n and have to be included into the 
trace over flavor space. 

Inserting Eqs. (58) and (62) into Eq. (57), we obtain for the general polarization tensor 



-Trfr^o^ A- It 70 7^ A£ 



+Tr 
+Tr 



KMln^rl-yoM^oKiK: 2 } <e 2 (ku k 2 ,Po) 
r^oA^To^A^r^n^J < ea (*i,fe,Po)} • (65) 



In the following we focus on the special case po = and p — > 0, i.e, k 2 — > ki = k. In this limit, the traces only depend 
on k and the functions v and w only on k = |k|. Thus, the c? 3 k integral factorizes into an angular and a radial part. 
With the abbreviations 



J2n) 3 



ao,eie2 



-^{Tr[r^ 70 nA k ei r^ 70 ^A k e2 ]+Tr T'^WPu f ll WPi-\ 



]}• 



+ Tr 



(66a) 



(66b) 



we can write the polarization tensor as 



IC(0) ^ limn^(0,p) = E E [ V abZ 2 l™/ *fc 2 ^ e2 (fc 1 ,fc 2 ,0) + WX ( : i S e2 lim| dk fc 2 < 62 fc 2 , 0) 



ei,e2 r,s 



(67) 

Note that only the angular integrals over the color, flavor, and Dirac traces, defined by V and W, depend on the 
symmetries of the order parameter and thus have to be calculated separately for each phase. Therefore, we first 
consider the dk integrals which are the same for all cases. 

In order to see how the two different quasiparticle excitations branches (labelled by r, s), as well as normal and 
anomalous propagation of the respective excitations (represented by the functions v, w) contribute in the final expres- 
sions for the polarization tensors in the zero-energy, zero-momentum limit, it is convenient to define the quantities 
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^ lim J dk k 2 [v r + \ (h , fc 2 , 0) + v r _ s _ (fci , fc 2 , 0)] , 
-llim [dkk 2 [v r + s _(ki,k 2 , 0) + v r _ s + (h,k 2 M , 
^ lim J dk k 2 [w r + s + (k u k 2 , 0) + w1 s _(fci, k 2 , 0)] , 
^lim [dkk 2 [w r + s _(k u k 2 ,0)+wL s + (k 1 ,k 2 ,0)] . 

p^O J 



(68a) 
(68b) 
(68c) 
(68d) 



These quantities are dimensionless since v™ e2 (fci , k 2 , 0) and w™ e2 (fci, k 2 , 0) have the dimension [1/energy] (cf. the 
definitions in Eqs. (59) and (63)). Combining particle-particle (ei = e 2 = +) and antiparticlc-antiparticlc (e\ — e 2 — 
— ), as well as particle-antiparticlc (ei = — e 2 = ±) contributions is possible since the corresponding integrals V, W 
multiplying these terms turn out to be the same. In the definitions of v rs , v rs and w rs , w rs we divided by the square 
of the quark chemical potential fi 2 in order to make these quantities independent of the quark flavor. This will be 
convenient for the results of the spin-one phases, where the formation of Cooper pairs is possible for different chemical 
potentials for each quark flavor. In App. A we present the calculation of the relevant integrals defined in Eqs. (68). 
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11 

V 


22 

D 


12 

V 


21 

V 


-11 
V 


-22 


-12 
V 


-21 


11 

w 


22 

W 


w 12 


W 21 


w 11 


-22 
W 


w 12 


w 21 


T = 


2SC, polar 


i 

2 


-1 




1 

2 


1 

2 


1 

2 











CFL, CSL 






1 

2 




1 

2 


1 

8 


1 

2 


iln2 





T>T C 


all phases 






1 




1 

2 






TABLE II: Leading order results for the integrals denned in Eqs. (68). The indices 1 and 2 correspond to the two gaps of 
each phase (the second one vanishing in the 2SC and the polar phase), while v corresponds to the trace over quasiparticle 
propagators G and w to the trace over anomalous propagators 3 . The fields with no entry indicate that these values do not 
occur in the calculations. 





M k 




vl 


2SC 




J| 


i-Ji 


CFL 


J 1 75 


|[(J-I) 2 -1)] 


±[4- (J- I) 2 ] 


polar 


J 3 [fc z +7l(k)] 


Ji 


i-Ji 


CSL 


J-[k + 7± (k)] 


i[k + 7± (k)][k- 7± (k)] 





TABLE III: Relevant color-flavor-Dirac matrices for the calculation of the Debye and Meissner masses in a given color- 
superconducting phase. The matrix .Mk reflects the symmetries of the various gap matrices. For the definition of the projectors 
P^' 2 see Eq. (f6). In color space, we use the matrices (Ji)jk = —icijk, (i,j,k = 1,2,3); in flavor space, we use {U)jk = —itijk 
and the second Pauli matrix T2. In Dirac space, we defined 7 ± (k) = 7 — k7 ■ k. 

As in Rcfs. [17, 18], we neglect the antiparticle gap, (f>~ ~ 0, and compute the integrals up to leading order. In Table 
II we collect all results. Some integrals vanish since they are proportional to a vanishing gap. This is the case for 
w rs = (for all r, s), since these integrals are proportional to at least one antiparticle gap. Clearly, for T > T Cl the 
gap vanishes in all phases and thus w rs = w rs = 0. 

In order to discuss the traces and the angular integral in Eq. (67), we have to distinguish between the several 
phases since the special form of the gap matrix, Eq. (13), is explicitly involved. Therefore, in the following sections 
we discuss the 2SC, CFL, polar, and CSL phases separately and compute the Debye and Meissner masses for photons 
and gluons in each phase. For the 2SC phase [17] and the CFL phase [18, 19], the results for the gluons are already 
known. Also the masses of the rotated gauge bosons, where the rotation is given by the new generator Q, cf. Eq. (51), 
were considered for these two phases [13]. Nevertheless, we will briefly discuss also these cases, since we first want to 
establish our notation and second, we will show that for the 2SC phase, there is actually no mixing between electric 
gluons and photons, i.e., the longitudinal mixing angle 6 2 is zero. Consequently, there is no rotated electric photon. 

The basic symmetry properties of all considered phases have been discussed in detail in Ref. [10]. We show the 
relevant color-flavor-Dirac matrices in Table III. 



B. The 2SC phase 

1. Gluon polarization tensor (a,b < S) 
Inserting the matrices given in the second line of Table III into Eqs. (66), we obtain 

= -9 2 Vr T * V °\ T Z,-e 2 + Tr[Tj V r Tff V s ] } , (69a) 

= 9 2 {Tr[T a J 3 V r If J 3 V s ] W£_ ea + Tr[Tj J 3 V r T b J 3 V s ] W\ , C2 } , (69b) 

where the traces only run over color space and where we defined 

T eT,e 2 = J T^WtoA^^oA-] , (70a) 

*Ce a = / ^MVVA^VA?] . (70b) 

Here, the traces only run over Dirac space. We used the fact that the projectors V r ' s = V^ s do not depend on the 
quark momentum k and that the color and Dirac traces factorize. Furthermore, the trivial flavor trace was already 



13 



performed, yielding a factor 2. The angular integrals over Dirac traces, Eqs. (70), are easily evaluated, 

1 



T 



oo 



^ 1 - ea = 2^(l + eie 2 ), 



n-0i yj-iO 7J®^ 7J^ 

ei,C2 ei,e2 ei,— e2 ei- — £2 

^', e2 = <,-e 2 = ^(l-^ ie2 ), 



(71a) 
(71b) 
(71c) 



where i, j = 1,2, 3. For the evaluation of the color traces note that J3V 1 — J 3 and J3P 2 = 0. We find that V and W, 
given in Eqs. (69), are diagonal in the adjoint color indices a and b. 

(a) [i = v = 0. With Eq. (71a) we obtain after performing the color traces 



2.gV 



X + X 1 



f (« 12 +« 21 ) 



1 11 , 1 22 _ 1 11 

< 6 V ~ 3 6 1 " 



for a = 1,2, 3, 



for a = 4 - 7 
for a = 8 . 



(b) fj, = 0,v = i and ^ = i,v = 0. Due to Eq. (71b), 

n°; b (o) = nj° (0) = . 

(c) fi = i,v = j. Due to Eq. (71c), the polarization tensor is diagonal in spatial indices i, j. We obtain 

i(w n +2V 11 )- i(w n +2w) n ) for a = 1,2,3, 

i[(w 12 +i; 21 ) + 2(i} 12 +i) 21 )] for a = 4- 7 

k i(w 11 +2w 11 ) + i(w 22 + 2w 22 ) + i( W 11 +2u; 11 ) for a = 8. 



(72) 



(73) 



(74) 



2. Mixed polarization tensor (a < 8, 6 = 9 and a = 9, fe < 8,) 

For a system in the color-superconducting 2SC phase where quarks with the electric charges qi and 52 form Cooper 
pairs, the electric charge generator, introduced in Eq. (4), is given by Q — diag(gi, g 2 )- Thus, we obtain 



v a~f,eie2 



\ eg( qi + q 2 ) {Tr[T a V r V s ] T^_ e2 + Tr[T a T V r V s ] 7^ } , 



(75a) 
(75b) 



WIZiZ, = 2 e 5(9i+92){Tr[T a J 3 ^J 3 7'1<% 2 +Tr[rJj 3 ^J37' s ]^ l!e2 } . 
It is not difficult to show that the polarization tensor is symmetric under the exchange of photon and gluon indices, 

n£(o) = n£(o) . (76) 

Using J| = J 3 , all color traces in Eq. (75) reduce to Tr[T a jf] = 5 aS /\/3. Therefore, we obtain for the various Dirac 
components of the tensor: 



(a) /U = v = 0. 



Il™(0) = ( qi + q2 ) 



v 11 — v 22 — w 11 



for a = 1 - 7 
for a = 8 . 



(77) 



(b) fi = 0, v = i and fj, = i,v = 0. 



n£(o) - n£(°) - . 



(78) 
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(c) n = i,v = j. 

for a = 1 - 7 , 
(v 11 + 2V 11 ) - (v 22 + 2v 22 ) + (w 11 + 2w n ) for a = 



n%(0) = 5 ij (q 1+q2 )^£{ (79) 



3V37T 



3. Photon polarization tensor (a — b — 9) 

In this case, the tensors V, W are 

VfiSe* = \ eHq! + 4) Tr[P r V s ] (T_^ i; _ e2 + 7^) , (80a) 
= e^Tr^^ 8 ] «% 2 +W^ liea ) . (80b) 

Here, we performed the flavor traces Tr[Q 2 ] = q\ + q 2 and Tr[QT 2 QT 2 ] = 2q\q 2 . After performing the color traces and 
the sums over e±, e 2 and r, s, the results for the different components are as follows. 



(a) n = v = 0. 



n°° 7 (o) = e ^4- M + qI) (2« n + v 22 ) 4q iq2W ^] . (si) 



(b) (i = 0,f = i and fj, = i,v = 0. 

n°;(o) = n;° 7 (o) = o . (82) 

(c) n = i,v = j. 

2 2 

n 77(°) = S H^2 {(«i 2 + fe 2 ) [2(« U + 2« n ) + (« 22 + 2v 22 )] + A qi q 2 (w 11 + 2W 11 )} . (83) 

C. The CFL phase 

1. Gluon polarization tensor (a,b < S) 
With the matrix Mk and the projectors V^' 2 for the CFL phase, given in Table III, Eqs. (66) become 

V„TeTe 2 = l 9 2 {Tr[T a r r T b r s }TZ.- e2 +TATlV r T^V s }T e % 2 } , (84a) 
W»Te" a = lg 2 {Tr[T a J-ir r T^J-ir s ]UZ- e2 +Tr[T^3-ir r T b 3-ir s }U^ eue2 } ■ (84b) 

As in the 2SC phase, the projectors V r , V s do not depend on momentum. Consequently, the angular integrals defined 
in Eq. (70) also appear in the CFL phase. Therefore, also in the CFL phase the (Oi) and (iO) components of the 
polarization tensor vanish, and the (ij) components are proportional to S lJ . Contrary to the 2SC phase, color and 
flavor traces cannot be performed separately, since the projectors V ', V s are nontrivial matrices both in color and in 
flavor space. In order to perform the color-flavor trace, one uses the relations Tr\T a 'P 1 T b T' 1 ] = and J • IV 1 — —2 V 1 . 
The polarization tensor is not only diagonal in color, but also has the same value for all eight gluons. One finally 
obtains the following expressions for the gluon polarization tensors. 

(a) fj, = v = 0. 

n°°6(0) - 5 ab ^ [(v 12 + v 21 ) + 7v 22 + 2(w 12 + w 21 ) + 2w 22 ] . (85) 
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(b) /j, = 0, v = i and [i = i,v = 0. 



n° a l(o) = Km = o . 



(86) 



(c) fj, = i,u = j. 



= S ll S ab ^£,[^ 2 +v 2 ^ + 2^ 2 +v 2 ^ + 7(v 22 + 2v 22 ) 



-2{w iz + w 2L ) - A(w L2 + w 2L ) - 2{w 22 + 2w 



(87) 



2. Mixed polarization tensor (a < 8,b = 9 and a = 9, b < 8) 



To compute the mixed polarization tensors, we need the electric charge generator Q. Since we consider a system 
with three quark flavors of electric charges qi, q 2 , (ft, we have Q = diag(gi, q 2 , (73). In the final result we will insert 
the charges for u, d, and s quarks. We obtain 



= \eg{TT[T a V r QV s ]T^ u _ e2 +Tr[TjV r QV s ]TZ 2 } 



1 



eg {Tr[T a J-IV r QJ-IV s }UZ-e 2 +^[TlJ-lV r QJ-IV s ] Wt^} ■ 



First we note that Eq. (76) also holds for the CFL phase. With the help of the relations 

Tr[T a V 1 QV 1 }=0 , Tr[T a Q] = Tr[TjTr[Q] = , 

and 

Tr[T a J ■ IQ J • I] = 3 Ti[T a V 1 Q] = 3Tr[T a V 1 Q3 ■ I] = S a3 ±( qi - q 2 ) + <*a8^=(gi + 92- 2q 3 ) 

we obtain the following results, 
(a) = v = 0. 



(88a) 
(88b) 

(89) 

(90) 



6ir 2 



(Qi - 92) [(v 12 + v 21 ) - 2v 22 + 2(w 12 + w 21 ) - 7w 22 )] 



for a =1,2,4- 7 
for a = 3 , 



k -±(qi + q 2 - 2q 3 ) [(v 12 + v 21 ) - 2v 22 + 2(w 12 + w 21 ) - 7w 22 )] for = 8. 
(b) (i = 0, v = i and fj, = i, v = 0. 



n«(o) = ir a ° 7 (o) = o 



(91) 
(92) 



(c) n = i,u = j. 



U%(0) = S* 



18tt 2 



(«i - 92) [(v 12 + v 21 ) + 2(v 12 + v 21 ) - 2v 22 - Av 22 
-2(w 12 + w 21 ) - A(w 12 + w 21 ) + 7w 22 + Uw 22 ] 

^(gi + 92 - 2 93 ) [(v 12 + v 21 ) + 2(v 12 + v 21 ) - 2v 22 - Av 
-2{w 12 + w 21 ) - A(w 12 + w 21 ) + 7w 22 + Uw 22 ] 



for a =1,2,4-7, 
for a = 3 , 
for a = 8 . 



(93) 
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3. Photon polarization tensor (a — b — 9) 

In this case, the tensors V, W read 

= \e 2 Tr[QV r QV s ] (T^^+TgJ , (94a) 
WfiZei = ^ 2 Tr[QJ.I7"-QJ-I^] «% 2 +UH v eu J . (94b) 
In order to perform the color-flavor traces, we abbreviate the following sums over the three quark charges, 

3 3 

x= IriQm ' V = X] q n ■ ( 95 ) 

n,m-l 7i—l 

Then, with Tr[Q V 1 Q V 1 } = x/9, Tr[Q 2 l c ] = 9Tr[QV 1 Q] = 3y (where l c is the unit matrix in color space), and 
Tr[Q J • IQ J • I] = 3Tr[Q J • \V X Q J • I] = 2{x - y), we obtain the following results. 

(a) n = v = 0. 

n °° 7 (0) = e ^-[ xv ^ + (3y-x)(v 12 +v 21 ) + (21y + x)v 22 



4xw LL + (6y-2x)(w LZ + w ZL ) + (6y - I0x)w z2 \ . (96) 



In a system of d, s, and u quarks we have the electric charges q\ = q 2 = —1/3 and = 2/3. In this case, 

x = , y=l. (97) 
Inserting these values into Eq. (96), the result becomes proportional to the gluon polarization tensor, Eq. (85), 

If! 

3? 



n°° (o) = -^nZ(o) ■ (98) 



(b) fj, = 0, v = i and ji = i, v = 0. 

(o) = n; 7 (o) = o . (99) 

(c) n = i,v = j. 

n«(0) = S tJ ^{x(v 11 +2v 11 ) + (3y-x)[v 12 +v 21 +2(v 12 +v 21 )] 
+ (21y + x)(v 22 + 2v 22 ) + Ax (w 11 + 2m 11 ) 

- (6y - 2x) [w 12 + w 21 + 2(w 12 + w 21 )] - (6y - I0x){w 22 + 2w 22 )} . (100) 

Again, using the quark charges of d, s, and u quarks that lead to Eq. (97), we obtain a result proportional to 
that given in Eq. (87), 

n 77(°) = |4 n aa(0) • (101) 

9 

D. The polar phase 

In the polar phase, Cooper pairs are formed by quarks of a single flavor and carry total spin one, J = 1. This phase 
is defined by the matrix Mk = Jz{k z +7^(k)], cf. Table III. The spin of the Cooper pair is aligned with the spatial z 
direction which means that rotational 5*0(3)7 is broken to 5*0(2)7, cf. Table I. As in the 2SC phase, the condensate 
also points in a fixed color direction. Physically, this means that quarks of one color, say blue, remain unpaired. This 
similarity to the 2SC phase can also be seen in the projectors V 1 ' 2 which are identical in both phases. Note that 
V 1 ' 2 in the polar phase do not depend on k although Mk does. In the following we consider a system of quarks with 
Nf different flavors where each quark flavor forms Cooper pairs separately. Each quark flavor has a separate electric 
charge, qi, . . . , qpf f , and chemical potential, . . . , Hn s ■ 
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1. Gluon polarization tensor (a, b < 8) 

In this case, we have 

VKa = lg^Tr[T a V r T b V s ]T^ u _ e2 +TT[Tjr r T^V s ]T^} , (102a) 

w%Z. = lg 2 {nT a j 3 r r T^j 3 r s ]u^_ e2 + Tr[r Q T j 3 v r nj 3 r s ]u^ e2 } , (102b) 

where 

UZe* =~f ^ Tr ^ (** - ^)) A? 7^ (fc z - 7 i(k)) A£»] . (103) 

Since 

Tr[y (fc* - 7 i(k)) A k x Y (k z - 7l(k)) A k fi2 ] = Tr[ 7 ^ A k x 7" A k £2 ] , (104) 

we find 

■ (105) 

Consequently, the gluon polarization tensor in the polar phase is almost identical to that in the 2SC phase. The only 
difference is the flavor trace which here yields a factor J2n=i ^n- Therefore, Eqs. (72) - (74) hold also for the polar 
phase after replacing 2^ 2 with Yl,nLi Ai- 

2. Mixed polarization tensor (a < 8,b = 9 and a — 9, b < 8) 

Obviously, also for the mixed polarization tensor, the Dirac and color part is identical to the 2SC case. Consequently, 
in Eqs. (77) - (79) one has to replace the factor (qi + (ft)^ 2 by Yln=i In^n in order to obtain the results for the polar 
phase. 

3. Photon polarization tensor (a = b = 9) 

Analogously, for the photon polarization tensor, the 2SC results given in Eqs. (81) - (83) are valid for the polar 
phase with the following modifications. Since in the 2SC phase there is a nontrivial flavor structure in the matrix 
Alk, there arc two different flavor factors in Eqs. (81) - (83), namely (q 2 + q 2 )n 2 and 2qiq2/J, 2 . Replacing each of 
these factors by the common factor Yln=i InA yields the corresponding results for the polar phase. 

E. The CSL phase 

As in Sec. IV D, we consider a system of Nj quark flavors, each flavor forming Cooper pairs separately. In the 
CSL phase, all quark colors acquire two gapped excitation energies. This is similar to the CFL phase, where there is 
also a quasiparticlc excitation branch with an energy gap <j> and a second one with an energy gap 2(f>. Formally, this 
structure has its origin in the two nonzero eigenvalues of the operator 7oA4 k .Mk7oj cf. Eqs. (15) and (16). In the 
case of the CSL phase, the matrix Aik as well as the projectors V^' 2 depend on the direction of the quark momentum 
k. Due to color-spin locking, the spatial components of the vector k are aligned with the three color directions red, 
green, and blue. 
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a 


V 00 ' 11 


-,,00,12 _ -,,00,21 


-, ,00,22 


1,3,4,6,8 





(1 + eie 2 )/(127r 2 ) 


(1 + eie 2 )/(127r i ) 


t 2,5,7 


(1 + eie 2 )/(187r ) 


(1 + eie 2 )/(367r ) 


5(1 + eie 2 )/(367r ) 




r v aa, ei&2 


Vy00,12 _ yyOO.21 


w oo, 22 


1,3,4,6,8 





(l + eie 2 )/(67r 2 ) 


r -(l + ei e 2 )/(127r 2 ) 


2,5,7 


2(l + eie 2 )/(97r 2 ) 


-(l + e ie2 )/(187r 2 ) 


5(l + eie 2 )/(367r 2 ) 



TABLE IV: (00) components of the tensors V, W for the gluon polarization tensor in the CSL phase, defined in Eqs. (106). 
All tensors are diagonal in color space. 



1. Gluon polarization tensor (a, b < S) 
Inserting the matrices Mu and P^' 2 from Table III into Eq. (66) yields 

KZU = y / ^{ Tr b"ioT a vi A k ei rionr^} + ( ei , 2 - -e ll2 , r Q , b - r a T fc ) } , (ioea) 

WE* - ~ / 7 0{Tr[ 7 ^T o J.(k- 7 Jk))nA^7,T- r J.(k- 7± (k))^A- 

+ ( e i,2 -» -ei,2, T 0) (, -> T a T b ) | . 



(106b) 



The trace over the 12 x 12 color-Dirac matrices is more complicated than in all previously discussed cases. For 
the (00) components, both V and W are diagonal in color space. The diagonal elements are collected in Table IV. 
They are divided into two parts, one corresponding to the symmetric Gcll-Mann matrices, a — 1,3,4,6,8, the other 
corresponding to the antisymmetric Gell-Mann matrices, a — 2, 5, 7. For the (iO) and (Oi) components, we find 
(omitting the indices r, s and ei, e-i) 



V, 



ab y ab ' y ab ' y ab 



o . 



(107) 



Therefore, also for the CSL phase, the (Oi) and (iO) components of the gluon polarization tensor vanish. 

For the (ij) components, we find that V and W are neither diagonal with respect to color a, b nor with respect 
to spatial indices. Nevertheless, after inserting the k integrals from Table II, the gluon polarization tensor becomes 
diagonal, i.e., n^ b (0) ~ S a bS l: > . The reason for the cancellation of all non-diagonal elements are the following properties 
of V and W, 



67T 2 



W a6,eie 2 ~ 1 _ eie2 for ^ 3\ a ; b arbitrary), and for (i = j; a ^ b) . 
Taking into account the trace over flavor space, which is the same as in the polar phase, we obtain 
(a) fj, = v = 0. 



(108a) 
(108b) 



N, 



n=l 



(b) fi = 0, v = i and = i,v = 0. 



3{v 12 + v 21 ) + Sv 22 
+6(w 12 + w 21 ) - 3w 22 

2V 11 + (v 12 +v 21 ) + 5v 22 
, +8W 11 - 2(w 12 + w 21 ) + 5w 22 



nS(o) = (o) = o . 



for a=l,3,4,6,8 



for a = 2, 5, 7 . 



(109) 



(110) 



(c) fj, = i, v = j. The general result (keeping all functions v, w) is a complicated 24 x 24 matrix and therefore 
not shown here. However, as stated above, the polarization tensor is diagonal after inserting the values for the 
functions v and w, and the result for the Debye and Meissner masses will be given in the next section. 
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2. Mixed polarization tensor (a < 8, b — 9 and a = 9, b < 8) 



Here we have 



' / a-y,eie 2 



f / ^{ Tr [^^ r ^kA k ei 7 ,y 7oQKA k e2 ]+K2-- e i, 2 ,T a ^rJ)}, (111a) 

+ (ei, 2 -» -ei, 2 , T a - Tj) } , (111b) 
where Q = diag(<?i, . . . , gjvj.)- All (00), (iO), and (Oi) components of these integrals vanish, 

yOO = yOi = ytO = yyOO = yyiO = yyOi = Q 

K a7 K a7 K a7 ' v 07 '"07 ' K a7 u • \ J 

(Here, we again omitted the indices r, s and e\, e 2 .) The (ij) components of V and W are nonvanishing and, as for 
the gluonic CSL case, we do not show them explicitly. However, the final result, i.e. the polarization tensor in the 
considered limit, vanishes because of the following relations, 

= 0, (113a) 

rs 

WjftU ~ l-eie 2 . (113b) 

(a) /x = i> = 0. According to Eq. (112), we have 

n°° 7 (o) = o . (in) 

(b) = 0, v = i and [i = i, v = 0. 

Il£(0) - H«(0) = . (115) 

(c) (i = i,v = j. This polarization tensor has a complicated structure in terms of v, w and is thus not given here. 
However, as stated above, the final result is zero. 



3. Photon polarization tensor (a — b — 9) 
For the photon polarization tensor in the CSL phase we need 

= y J -^{Tr [/ToQ^Ar 1 / 7oQ W 2 ] + («i ; 2 - ^1,2) } , (116a) 

W#£ ea = ~ |^{Tr[ 7 ^QJ.(k- 7± (k))KA k 1 7 ,QJ-(k-7 i (k))nA k e ^ 

+ (ei, 2 - -ei, 2 ) } ■ (H6b) 

The results for the (00) and (ij) components are given in Tables V and VI. The (Oi) and (iO) components vanish (for 
all r, s, ei, e 2 ), 

v;° 7 = v°; = w;° 7 = w°; = o . (117) 

One obtains the following photon polarization tensor, 
(a) fj, = v = 0. 

E£(0) = £ « 4(« U + 2- 22 - - 2 W 22 ) . (118) 



20 



v oo,n 

1/ 77, e l e 2 


y00,12 _ -1)00,21 
1/ 77,eie 2 1/ 77,ei e 2 


■1)00,22 
^77,61 e 2 


(l + eie 2 )/(27r 2 ) 





(l + eie 2 )/7r 2 


vyoo.n 

77! e l e 2 


vy00,12 _ yuOO.21 
r K 77,eie2 ^ 77> e l e 2 


yyOO.22 


-2(1 + eie 2 )/7r 2 





-(l + eie 2 )/7r 2 



TABLE V: (00) components of the tensors V, W for the photon polarization tensor in the CSL phase, defined in Eqs. (116). 



"77,eie 2 


vw. 12 — V''' 21 

^77,eie 2 K 77, e l e 2 


■\jij,22 

K 77,eie 2 


(5' 3 (ll + 7eie 2 )/(547r 2 ) 


8<5 lJ (l-eie 2 )/(277r 2 ) 


(5 I3 (19-eie 2 )/(277r 2 ) 


' ^77,eie 2 


yyij, 12 _ yyij, 21 
' *77,eie 2 ' *77,eie 2 


yyij,22 

' ^77,eie 2 


2<f J (7 + lle ie2 )/(277r 2 ) 


16(5 lJ (l-eie 2 )/(277r 2 ) 


-£"(1 - 19eie 2 )/(277r 2 ) 



TABLE VI: (ij) components of the tensors V, W for the photon polarization tensor in the CSL phase, defined in Eqs. (116). 

(b) [/, = Q,i' = i and fi = i, v = 0. 

n£(o) = n; 7 (o) = o . (119) 

(c) [l = i,V = j. 

II«(0) - ^^ g X^[^ 11 + 2i) 11 + 16(S 12 + z; 21 ) + 18 W 22 + 20i; 22 

n=l 

+ 36w n -8w n +32(w 12 +w 21 ) + 18w 22 - 20u; 22 ] . (120) 

V. MIXING AND SCREENING 

In this section we use the results of the previous section to calculate the Debye and Meissner masses. We insert 
the numbers from Table II into the results for the polarization tensors n°°(0) and n^ b (0) given in Sees. IV B, IV C, 
IV D, IV E and use the definitions of the screening masses given in Sec. II A (cf. Eqs. (39) and comment below this 
equation). 

We distinguish between the normal-conducting and the superconducting phase. The screening properties of the 
normal-conducting phase are obtained with the numbers given in Table II for temperatures larger than the critical 
temperature for the superconducting phase transition, T > T c . For all a,b < 9, they lead to a vanishing Meissner 
mass, i.e., as expected, there is no Meissner effect in the normal-conducting state. However, there is electric screening 
for temperatures larger than T c . Here, the Debye mass solely depends on the number of quark flavors and their electric 
charge. We find (with a < 8) 

T>T C : ml aa = 3N f 9 -^- , = , m 2 , ^ = 18 £ g 2 . (121) 

n 

Consequently, the 9x9 Debye mass matrix is already diagonal. Electric gluons and electric photons are screened. 

The masses in the superconducting phases are more interesting. The results for all phases are collected in Table VII 
(Debye masses) and Table VIII (Meissner masses). The physically relevant, or "rotated", masses are obtained after a 
diagonalization of the 9x9 mass matrices. We see from Tables VII and VIII that the special situation discussed in 
Sec. II B applies to all considered phases; namely, all off-diagonal gluon masses, a, b < 8, as well as all mixed masses 
for a < 7 vanish. Furthermore, in all cases where the mass matrix is not diagonal, we find mg 7 = m^m^. Therefore, 
the rotated masses frigs and m 77 (which are the eigenvalues of the mass matrix) are determined by Eqs. (45). The 
electric and magnetic mixing angles 9d = 2 and 9m = 9\ are obtained with the help of Eqs. (44). Remember that 
the indices 1 and 2 originate from the spatially transverse and longitudinal projectors defined in Eqs. (18). They were 
associated with the Meissner and Debye masses in Eqs. (39). We collect the rotated masses and mixing angles for all 
phases in Table IX. 

Let us first discuss the spin-zero cases, 2SC and CFL. In the 2SC phase, the gluon masses are obtained from Eqs. 
(72) and (74). Due to a cancellation of the normal and anomalous parts, represented by v and w, the Debye and 
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m D,ao 


2 2 
m D,ai = m D,~/a 
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a 


1 2 3 4 5 6 7 8 


1-7 8 
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2SC 












2e 2 
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3<ff 2 
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3 5 2 





18<? 2 e 2 
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TABLE VII: Zero-temperature Debye masses. All masses are given in units of N f p? j '(67r 2 ) , where Nf = 2 in the 2SC phase, 
Nf = 3 in the CFL phase, and Nf = 1 in the polar and CSL phases. We use the abbreviations £ = (21 — 81n2)/54, 
a = (3 + 41n 2)/27, and (3 = (6 - 4 In 2) /9. 
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TABLE VIII: Zero-temperature Meissner masses. All results are given in the same units as the Debye masses in Table VII. 
The abbreviations of Table VII are used. 

Meissner masses for the gluons 1,2, and 3 vanish. Physically, this is easy to understand. Since the condensate picks 
one color direction, all quarks of the third color, say blue, remain unpaired. The first three gluons only interact with 
red and green quarks and thus acquire neither a Debye nor a Meissner mass. We recover the results of Rcf. [17]. For 
the mixed and photon masses we inserted the electric charges for u and d quarks, i.e., in Eqs. (77), (79), (81), and 
(83) we set qi = 2/3 and q2 = —1/3. Here we find the remarkable result that the mixing angle for the Debye masses 
is different from that for the Meissner masses, / 6m- The Meissner mass matrix is not diagonal. By a rotation 
with the angle 9m, given in Table IX, we diagonalize this matrix and find a vanishing mass for the new photon. 
Consequently, there is no electromagnetic Meissner effect in this case. This fact is well-known [12]. The Debye mass 
matrix, however, is diagonal. The off-diagonal elements rr? D % vanish, since the contribution of the ungapped modes, 
corresponding to the blue quarks, v 22 , cancels the one of the gapped modes, w 11 — w 11 , cf. Eq. (77). Consequently, the 
mixing angle is zero, Od = 0. Physically, this means that not only the color-electric eighth gluon but also the electric 
photon is screened. Had we considered only the gapped quarks, i.e., v 22 = in Eqs. (72), (77), and (81), we would 
have found the same mixing angle as for the Meissner masses and a vanishing Debye mass for the new photon. This 
mixing angle is the same as predicted from simple group-theoretical arguments, Eq. (55). The photon Debye mass in 
the superconducting 2SC phase differs from that of the normal phase, Eq. (121), which, for q\ = 2/3 and 92 = — 1/3 
is m A 77 = 5A7eV/(67r 2 ). 

In the CFL phase, all eight gluon Debye and Meissner masses are equal. This reflects the symmetry of the condensate 
where there is no preferred color direction. For the mixed and photon masses, we used Eq. (97), i.e., we inserted 
the electric charges for u, d, and s quarks into the more general expressions given in Sec. IV C. The results in 
Tables VII and VIII show that both Debye and Meissner mass matrices have nonzero off-diagonal elements, namely 
TOg 7 = m 2 s . Diagonalization yields a zero eigenvalue in both cases. This means that neither electric nor magnetic 
(rotated) photons are screened. Or, in other words, there is a charge with respect to which the Cooper pairs are 
neutral. Especially, there is no electromagnetic Meissner effect in the CFL phase, either. Note that the CFL phase is 
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TABLE IX: Zero-temperature rotated Debye and Meissner masses in units of Nf fj, 2 / (6-k 2 ) and mixing angles for electric and 
magnetic gauge bosons. The constants £, a, and /3 are defined as in Tables VII and VIII. 
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the only one considered in this paper in which electric photons are not screened. Unlike the 2SC phase, both electric 
and magnetic gauge fields are rotated with the same mixing angle 9 D = 6 M . This angle is well-known [12, 13]. Note 
that, due to the spectrum of the matrix 7o-Mj c .Mk7o = (J ■ I) 2 , there are two gapped branches. Unlike the 2SC phase, 
there is no ungapped quasiparticle excitation branch. This is the reason why both angles 9d,&m coincide with the 
one predicted in Eq. (55). 

Let us now discuss the spin-one phases, i.e., the polar and CSL phases. For the sake of simplicity, all results in 
Tables VII, VIII, and IX refer to a single quark system, Nf — 1, where the quarks carry the electric charge q. After 
discussing this most simple case, we will comment on the situation where Nf > 1 quark flavors separately form Cooper 
pairs. The results for the gluon masses show that, up to a factor Nf, there is no difference between the polar phase 
and the 2SC phase regarding screening of color fields. This was expected since also in the polar phase the blue quarks 
remain unpaired. Consequently, the gluons with adjoint color index a = 1,2,3 are not screened. Note that, due to 
Eq. (105), the spatial z-direction picked by the spin of the Cooper pairs has no effect on the screening masses. As 
in the 2SC phase, electric gluons do not mix with the photon. There is electromagnetic Debye screening, which, in 
this case, yields the same photon Debye mass as in the normal phase, cf. Eq. (121). The Meissner mass matrix is 
diagonalized by an orthogonal transformation defined by the mixing angle which equals the one in Eq. (55). 

In the CSL phase, we find a special pattern of the gluon Debye and Meissner masses. In both cases, there is a 
difference between the gluons corresponding to the symmetric Gell-Mann matrices with a = 1,3,4,6,8 and the ones 
corresponding to the antisymmetric matrices, a = 2,5, 7. The reason for this is, of course, the residual symmetry 
group SO(3) c+ j that describes joint rotations in color and real space and which is generated by a combination of the 
generators of the spin group 5*0(3)7 an d the antisymmetric Gell-Mann matrices, T%, T5, and T7. The remarkable 
property of the CSL phase is that both Debye and Meissner mass matrices are diagonal. In the case of the Debye 
masses, the mixed entries of the matrix, rri£,^, are zero because of the vanishing traces, Eq. (112), indicating that 
pure symmetry reasons are responsible for this fact (remember that, in the 2SC phase, the reason for the same fact 
was a cancellation of the terms originating from the gapped and ungapped excitation branches) . There is a nonzero 
photon Debye mass which is identical to that of the polar and the normal phase which shows that electric photons are 
screened in the CSL phase. Moreover, and only in this phase, also magnetic photons are screened. This means that 
there is an electromagnetic Meissner effect. Consequently, there is no charge, neither electric charge, nor color charge, 
nor any combination of these charges, with respect to which the Cooper pairs are neutral. This was also shown in 
Sec. Ill where we argued that in the CSL phase there is no nontrivial residual U(l) em , cf. Table I. This feature of the 
CSL phase was already discussed in Ref. [20]. 

Finally, let us discuss the more complicated situation of a many- flavor system, Nf > 1, which is in a superconducting 
state with spin-one Cooper pairs. In both polar and CSL phases, this extension of the system modifies the results 
in Tables VII, VIII, and IX. We have to include several different electric quark charges, q\, ■ ■ ■ ,qN f , and chemical 
potentials, . . . ,/ijVf , in a way explained below Eq. (105) and shown in the explicit results of the CSL phase, 
Sec. IV E. In the CSL phase, these modifications will change the numerical values of all masses, but the qualitative 
conclusions, namely that there is no mixing and electric as well as magnetic screening, remain unchanged. In the 
case of the polar phase, however, a many-flavor system might change the conclusions concerning the Meissner masses. 
While in the one-flavor case, diagonalization of the Meissner mass matrix leads to a vanishing photon Meissner mass, 
this is no longer true in the general case with arbitrary Nf. There is only a zero eigenvalue if the determinant of 
the matrix vanishes, i.e., if m 2 M Sy — mM,s8 m M^f Generalizing the results from Table VIII, this condition can be 
written as 



Consequently, in general and for fixed charges q n , there is a hypersurface in the V/-dimensional space spanned by 
the quark chemical potentials on which there is a vanishing eigenvalue of the Meissner mass matrix and thus no 
electromagnetic Meissner effect. All remaining points in this space correspond to a situation where the (new) photon 
Meissner mass is nonzero (although, of course, there might be a mixing of the eighth gluon and the photon). Eq. 
(122) is trivially fulfilled when the electric charges of all quarks are equal. Then we have no electromagnetic Meissner 
effect in the polar phase which is plausible since, regarding electromagnetism, this situation is similar to the one-flavor 
case. For specific values of the electric charges we find very simple conditions for the chemical potentials. In a two 
flavor system with qi = 2/3, q2 = —1/3, Eq. (122) reads 




(122) 



m,n 



= 



(123) 



while, in a three-flavor system with qi = —1/3, q2 = —1/3, and q 3 = 2/3, we have 



(124) 
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Consequently, these systems always, i.e., for all combinations of the chemical potentials /i„, exhibit the electromagnetic 
Meissner effect in the polar phase except when they reduce to the above discussed simpler cases (same electric charge 
of all quarks or a one-flavor system). 

VI. SUMMARY AND CONCLUSIONS 

We have presented a calculation of the gluon and photon polarization tensors U.^(P), a, b = 1, . . . , 8, 7, in a color 
superconductor which yield, in the zero-energy, low-momentum limit, the Debye and Meissner masses for the gauge 
bosons. These masses indicate in which cases electric and magnetic fields are screened. Here, magnetic screening is 
equivalent to the Meissner effect. We have derived the gauge field propagators in a general treatment starting from 
the QCD partition function, ft has been shown that, in general, longitudinal and transverse modes of the gauge fields 
are mixed in a nontrivial way. An "unmixing" transformation leads to a modified longitudinal polarization tensor. 
However, in the static limit, this mixing is absent and longitudinal and transverse modes are well separated. 

We have shown that a nondiagonal polarization tensor causes a mixing of the gauge fields. In general, the tensor has 
to be diagonalized via a transformation generated by an orthogonal operator 0{P). Consequently, also the gluon and 
photon fields, A®, are transformed. The operator O(P) defines new gauge fields, A a , which are linear combinations 
of the original ones. This is in complete analogy to the standard model of electroweak interactions, where, due to 
spontaneous breaking of the electroweak symmetry SU(2) x U(l), the new fields are the gauge bosons of the weak 
and electromagnetic interaction. In this case, the three fields corresponding to the weak interaction attain a mass via 
the Higgs mechanism. Similarly, also in the case of the color-superconducting phase transition, some of the originally 
nine gauge fields can become massive (besides a possible mixing). In the framework of the formalism presented here, 
the masses of the new fields are the eigenvalues of the polarization tensor n^ fc (0), where the index i = 1, 2 corresponds 
to the transverse and longitudinal modes (Meissner and Debye masses, respectively). Thus, the number of vanishing 
eigenvalues of n^ fc (0) indicate the dimension of the unbroken subgroup of the color-electromagnetic product group 
SU(3) C x U(l) em . 

We have presented a general formalism to calculate the polarization tensor for different color-superconducting 
phases. We have explicitly computed its zero-energy, low-momentum limit for the 2SC, CFL, polar, and CSL phases. 
Parts of the results were already known in the literature, namely the gluon Debye and Meissner masses for the spin-0 
phases 2SC and CFL [17-19]. Our result for the photon Debye mass in the 2SC phase differs from that of Ref. [13] 
since our calculation shows that the photon-gluon mass matrix is already diagonal and thus electric gluons do not 
mix with the photon (while in Ref. [13] the same diagonal matrix was rotated and a "mixed mass" was obtained). 
The masses for the spin-1 phases have never been computed before. For the polar phase, we have shown that there 
is mixing between the magnetic gauge bosons but, as in the 2SC phase, no mixing of the electric gauge bosons. In a 
system of one quark flavor, this mixing leads to a vanishing Meissner mass. However, for more than one quark flavor, 
we have shown that, if the electric charges of the quark flavors are not identical, there is an electromagnetic Meissner 
effect in the polar phase, contrary to both considered spin-0 phases. For the CSL phase, we have found the remarkable 
result that, for any number of flavors, neither electric nor magnetic gauge fields are mixed. Since there is no vanishing 
eigenvalue of n* 6 (0), all eight gluons and the photon (electric as well as magnetic modes) become massive and there 
is an electromagnetic Meissner effect. This might affect the electromagnetic properties of a neutron star that has a 
core of color-superconducting quark matter in the CSL phase. We have already discussed these possible astrophysical 
implications in Ref. [20]. There we also argued that, in spite of a suppression of the gap of three orders of magnitude 
compared to the spin-0 gaps [9, 10], spin-1 gaps might be preferred in a charge-neutral system. The reason for that 
is that a mismatch of the Fermi surfaces of different quark flavors has no effect on the spin-1 phases, where quarks of 
the same flavor form Cooper pairs. 
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APPENDIX A: INTEGRALS OVER QUARK MOMENTUM k 

In this appendix, we prove the results shown in Table II, i.e., we calculate the integrals over quark momentum 
defined in Eqs. (68). In this calculation we will neglect the antiparticle gap, 0~ ~ 0. 
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1. Contributions from normal propagators at T — 



We start from Eqs. (59) which, for p = reduce to (cf. Eq. (61)) 
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-(1 - jVi, r - N 2>s ) + 



ei,r-e2, s + 66 



2ei,r-e2, s (ei,r + £ 2,s) v 2ei jT -e2, s (ei,r- — £2,s) 

Here, we inserted the definition of n i!r given in Eq. (60) and defined 

In the limit p — > 0, or, cquivalcntly, &i — > k 2 , we obtain for r = 1,2 



(JVl.r ~ ^2, s ) 
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2 / rffcfc' 



4e3 



■(1 - 2JV r ) 
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2e? de r 



(Al) 



(A2) 



(A3) 



where we used (iVi, r — ^2,r)/(^i,r — C2,r) — dN r /de r and £ = fc — /x. All quantities now depend on fci = fc and 
correspond to positive energies, ei = e2 = 1, since, due to <j>~ ~ 0, the terms corresponding to negative energies 
vanish. Therefore, we omit the index i = 1,2. Note that, up to sublcading order, the gap function does not depend 
on the index r, <pi ~ <p 2 = 4> [10]. We have to distinguish between the cases where A r ^ and where A r = 0. When 
A r 7^ 0, we obtain for zero temperature T = 0, where N r = 0, 



1 Z" 00 \ J,2 1 ft 1 

A 4 JO 4e r z JO 



3/2 1 

2 



(A4) 



Here, we restricted the fc-integration to the region < k < 2/j, since the gap function is strongly peaked around the 
Fermi surface, k — ji. In this region, we assume that the gap function does not depend on the momentum k. Since 
the eigenvalue A r cancels out, this result does not depend on r. Therefore it gives rise to the value of v 11 for all 
four considered phases and for v 22 in the cases of the CFL and CSL phases. When A r = 0, we obtain (for arbitrary 
temperature T) 



1 f c 

v rr = - 

M Jo 



dk k 



2 dNp 
~dk~ 



(A5) 



where Np = l/[exp(£_/T) + 1] is the Fermi distribution. Using fi » T, and the substitution ( = £/2T, this can be 
transformed to 



v rr ~-f 
Jo 



<K 



cosh £ 



-1 . 



(A6) 



Let us now compute v rs for r ^ s. Using Eq. (Al), we obtain for two nonvanishing eigenvalues Ai,A2 and at zero 
temperature 



- f 

M 2 io 



dk k 1 



£l£2 



2ei£2(ei + e 2 ) 



(A 2 - Ai)</> 2 



4+e 4+e 



£2 



ei 



(A7) 



Again, we neglected the integral over the region k > 2/x. In the last integral, the leading order contributions cancel. 
The subleading terms give rise to (A 2 — \i)(f> 2 /2. Therefore, the eigenvalues cancel and the result does not depend on 
Ai, A 2 ^ 0. Eq. (A7) holds for the CFL and CSL phases. 

When one of the excitation branches is ungapped (2SC and polar phase), for instance, A 2 = 0, we have 



v 12 =v 21 = - 



1 

A* 2 Jo 



dk k 1 



For T = 0, this becomes with the standard approximations 
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(A9) 
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where the last integral has already been computed in Ref. [17]. 

Next we compute v rs for the various cases. Setting the thermal distribution function for antiparticles to zero, we 
obtain with the definition in Eq. (68) 
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l^Jo 



dkk 2 



(l-n r )(l-N r ) 



N r n r _ + (1 - n 8 )(l - N s ) 



N.n. 



et - e 



€ s € r 



(A10) 



with ef = e^ r . Again, we first consider the situation where r = s. In this case, for T = and A r ^ 0, (defining 

6 = 6+) 



II 2 I 

h 1 Jo 



II 2 I 

M Jo 



dkk 2 
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dk k 



n r 



1 

2k 



e r + fc + M 
/u(e r + jj, — k) + X r <p 2 1 
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(All) 



Here, — l/(2fc) is a vacuum subtraction. The last integral has been computed in Ref. [17]. Note that the result does 
not depend on the value of A r ^ 0. Thus, it is valid for w 11 in the 2SC and polar phases as well as for both v 11 and 
v 22 in the CFL and CSL phases. 

For A r = 0, we obtain from Eq. (A10) at zero temperature and with the vacuum subtraction l/(2fc) 



a 2 i 
M Jo 



dk k 1 



@(k-fj,) 1 
2k 2k 



(A12) 



This result holds for v 22 in the 2SC and polar phases. 

Next, we discuss the case r ^ s. From Eqs. (A10) and (All) it is obvious, since the result for v rr did not depend 
on A r , that we get the same result for v rs with two nonvanishing eigenvalues Ai, A 2 ^ 0. Thus, in this case, 



v 12 = v 21 = 



(A13) 



For A2 = 0, we obtain at zero temperature 



v 12 = v 21 = 
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2k 
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(A14) 



where identical integrals as in Eqs. (All) and (A12) were performed. Consequently, also for the 2SC and polar phases, 
Eq. (A13) holds. 



2. Contributions from anomalous propagators at T = 



We start from Eq. (63) which, for p = 0, is 



w rs = 

eie 2 



2ei l7 -e2,s(ei,r + £2,s) 



(1 - JVi, r - N 2tS ) + 
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(N hr - N 2 , s ) 
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Obviously, since <f>~ ~ 0, we have for all r, s and all phases w rs — (for the definition of w rs cf. Eq. (68)). First, 
we calculate w rs for r — s. Neglecting the antiparticle gap and taking the limit k\ — > fc 2 , we obtain, using the same 
notation as above, e r = e\ T , 



n 2 1 

M Jo 



— 7t I dkk 2 



f-(l-2iV r ) 



<j> 2 dN r 



For A r 7^ and at zero temperature, where N r — 0, this expression reduces to 



7 2 S? kk %-\fo 5 ^ e+x ^ /2 -k- 
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(A17) 



Unlike in v rr , the eigenvalue A r does not cancel out and we obtain different results for A r = 4 and A r = 1. In the cases 
of the 2SC and polar phases we obtain w 11 — 1/2 whereas for the CFL and CSL phases, w 11 = 1/8 and w 22 = 1/2. In 
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the former two cases, the quantity w 22 does not occur in our calculation. Thus we can turn to the case where r ^ s. 
Here, we conclude from Eq. (A15) 
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For two nonvanishing eigenvalues Ai, A 2 , this reads at T = 
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For Ai = 4 and A2 = 1 (CFL and CSL phases) we obtain 



(A18) 



(A19) 



w 12 = w 21 = -ln2 . 
3 



(A20) 



The corresponding expression for the case where the second eigenvalue vanishes, A 2 = 0, docs not occur in our 
calculation. 



3. Integrals in the normal phase, T > T c 



For temperatures larger than the transition temperature T c but still much smaller than the chemical potential, 
T <C n, all integrals defined in Eqs. (68) are easily computed. Since there is no gap in this case, (f> — 0, all 
contributions from the anomalous propagators vanish trivially, w rs — w rs = 0. From Eq. (AI), which hold for all 
temperatures, we find with = 



11 22 12 
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and (with the vacuum subtraction 1/fc) 
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